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Abstract 

We study a hierarchy of integrable partial differential equations 
in 2 + 1 dimensions arising from the commutation of 2 - dimen- 
sional vector fields, and we construct the formal solution of the as- 
sociated Cauchy problems using the Inverse Scattering Transform for 
one-parameter families of vector fields recently introduced in [TJ [2]. 
Due to the ring property of the space of eigenfunctions, the inverse 
problem can be formulated in three distinguished ways; in particular, 
one formulation corresponds to a linear integral equation for a Jost 
eigenf unction, and another formulation is a scalar nonlinear Riemann 
problem for suitable analytic eigenfunctions. 



1 Introduction 

In this paper we study a hierarchy of integrable partial differential equations 
(PDEs) in 2 + 1 dimensions arising from the commutation of 2 - dimensional 
vector fields, and we construct the formal solution of the associated Cauchy 
problems using the Inverse Scattering Transform (1ST) for vector fields re- 
cently introduced in [lj [2] . 

The first nontrivial member of such hierarchy is the equation 

V X t + Vyy = VyV xx - V x V X y, V = V(x, ]J, t) G R, X , JJ , t G R, (l) 

introduced in [3] and studied in [3j HJ E] from different points of view, which 
admits the operator representation 

[L, M] = (2) 
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in terms of the 2- dimensional vector fields [5] 



L = d y + (p + v x )d x , , , 

M = d t + (p 2 +pv x -v y )d x , 

where p is a constant parameter. 

Equation (pQ) is the u = reduction of the recently introduced [T] inte- 
grable system: 

Uxt + Uyy = -{UU X ) X - V X U X y + VyU XX , U, V E R, X,l/,t E R, 

Vxt "T Vyy V>V XX ^x^xy ~t~ ^y^xx 

admitting the operator representation [L, M] = in terms of the following 
3-dimensional vector fields 

L^=d y + {p + v x )d x - u x d p , , . 

M = d t + (p 2 + pv x +u- v y )d x + {-pu x + u y )dp. 

We remark that the v — reduction of the system is the celebrated dis- 
persionless Kadomtsev-Petviashvili (dKP) (or Khokhlov-Zabolotskaya) equa- 
tion: 

u tx + Uyy + (uu x ) x = 0, u — u(x, y, t) E R, x,y,tER. (6) 

arising in various problems of Mathematical Physics, and intensively studied 
in the recent literature (see, f.i., |S1 H El El HDl H). 

The inverse scattering problem for one - parameter families of multi- 
dimensional vector fields has been recently introduced in p], and used to 
construct the formal solution of the Cauchy problem for a wide class of non- 
linear PDEs in arbitrary dimensions [T], including the heavenly equation of 
Plebanski [T2] , and for the system (Jl]) [2j . The Hamiltonian constraint on the 
associated spectral data leading to the heavenly equation and to the dKP 
reduction of (j3J) was also found in pQ and [2]. We remark that, in the case of 
Hamiltonian vector fields, an elegant and alternative integration scheme was 
already known in the literature [?]. 

We also remark that, since the hierarchy of PDEs investigated in this 
paper arises from the commutation of 2 - dimensional vector fields, according 
to the theory developed in [H 12], the associated inverse problem can be 
formulated in terms of a scalar nonlinear Riemann problem. 
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2 Hierarchy of PDEs associated with L 

Equation ([TJ) is the first nontrivial member (corresponding to n = 2) of the 
following hierarchy of PDEs in 2 + 1 dimensions arising from the commutation 
of 2 - dimensional vector fields: 



v tn + Q n v x = 0, n e N, 
where Q is the recursion operator 

Q = d~ l (v xx -d y - v x d x ) =v x - d~ l d y - 2d~ 1 v x d x . 
The first few flows are: 



Vto +v x =0, 
v tl -v y = 0, 

V X t2 Vyy VyV xx -\~ V x V X y 0, 

V XXt3 - (Vy + Vl) yy + V X y (Vy - Vl) + V XX d~ l (Vy + ^), 

This hierarchy admits the operator representation 



(7) 



(8) 



(9) 



L,M n 



0, nGN, 



(10) 



where L is defined in fl3]) and the 2-dimensional vector fields M n are defined 
by 



n-l 



fc=0 



(11) 



with: 



4 n +i = Q4 nJ > !<*<"-!, A™=v x . (12) 

To show it, we replace and ffTTl) into ffTU]) . obtaining the following poly- 
nomial (in p) equation: 



n-l 



n-l 



n-l 



V xt + V xx p» + ^/iW = J2p k ( A n-k)v + (P + ^) T,P k ( A n-k)*> ( 13 ) 



n— fc 

fc=0 / fc=0 



fc=0 



which must be satisfied Vp. Equating to zero all powers of p, we obtain 
equations (1T21) and (J7J). 
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Since Q does not depend explicitly on x, its Lie derivative along v x , the 
first vector field of the hierarchy, is zero; in addition, it is also possible to 
verify that Q is a Nijenhuis (hereditary) operator. Therefore all the flows ([7j) 
commute with each other [131 E] ■ 

We remark that the algebraic theory for the 2 - dimensional recursion 
operator Q is conceptually more similar to the theory for the 1 - dimensional 
recursion operators associated with systems like the Korteweg - de Vries and 
nonlinear Shrodinger equations [13], rather than to the theory for the 2 - 
dimensional recursion operators associated with systems like the Kadomtsev 
- Petviashvili and Davey Stewartson equations [15]. No results are known 
though, at the moment, on the Hamiltonian and bi-Hamiltonian [T5] charac- 
ter of the hierarchy (J7J). 



3 Inverse scattering transform 

Now we consider the Cauchy problem for the hierarchy (J7|) and, in particular, 
for equation ([I]), within the class of rapidly decreasing real potentials v. 

v(x,y,t n )^0, (x 2 + y 2 )^oo, uGR, (x,y)eR 2 , t n >0, (14) 

interpreting t n as time and the other two variables x, y as space variables. 
To solve such a Cauchy problem by the 1ST method developed in [U we 
construct the 1ST for the operator L, within the class of rapidly decreasing 
real potentials v, interpreting the operators M n as time operators. 

The operator representation ({TO]) implies the existence of common eigen- 
functions f(x,y,t,p) of L,M n (the Lax pair): 

Lf = 0, M n f = 0, neN. (15) 

Since the Lax pair (fT5l) is made of vector fields, the space of eigenfunctions 
is a ring: if fx, f'2 are two solutions of the Lax pair (1151) . then an arbitrary 
different iable function F(f\, f 2 ) of them is a solution of ffTBI) . 



3.1 Direct problem 

The localization (|I4"I) of the potential v implies that, if / is a solution of 
Lf = 0, then 

f(x,y,p) -> f±(£,p), y-^±oc, « 16 . 
f := x-py; 
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i.e., asymptotically, / is an arbitrary function of £ = x — py and p. 

A central role in the theory is played by the real Jost eigenfunction 
<p(x, y,p), the solution of Lip = uniquely defined by the asymptotics 

<p(x,y,p)^£, y -> -oo, (17) 

and equivalently characterized by the linear integral equation ip = £ + 
G(— v x ip x ), for the Green's function G(x,y,p) = 9(y)5(x — py). 

The y = +00 limit of ip defines the natural scattering datum a for L: 



v- 



lim tp{x,y,p) = S(£,p)=£ + <r(£,p). (18) 



The direct problem is the transformation from the real potential v, func- 
tion of the two real variables (x,y), to the real scattering datum a, function 
of the two real variables (£,p). Therefore the mapping is consistent. In the 
small field limit, this mapping reduces to the Radon transform 



°~(£,P) = - / vz(€ + py,y)dy. (19) 
Jr 

A crucial role in the 1ST theory for the vector field L is also played by 
the analytic eigenfunctions il>±(x, y,p), the solution of Lif)± = satisfying 
the integral equations 

ij>±(x,y,p) = -J* R 2 dx'dy'G ± (x - x',y - y' ,p)v x <(x' ,y')ip± x ,(x' ,y' ,p) + f , 

(20) 

where G± are the analytic Green's functions 

G ± (x,y,p) = ±— ] (21) 

The analyticity properties of G±(x, y,p) in the complex p - plane imply that 
ip + (x,y,p) and ip-(x,y,p) are analytic, respectively, in the upper and lower 
halves of the p - plane, with the following asymptotics, for large p: 

M^y, P )=^-^^- + o(\), \ P \»i. (22) 

V \P J 

It is important to remark that the analytic Green's functions (1211) exhibit 
the following asymptotics for y — ► ±00: 

G±(x -x',y- y',p) -> ± — L^ , y -> +00, 
G ± (x - x', y - y 7 , p) -> ± 2 ^ K _ 1 g , ±ie] , y -> -00, 
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entailing that they = +00 asymptotics ofip + andijj_ are analytic respectively 
in the lower and upper halves of the complex plane £, while the y = —00 
asymptotics of ip + and ip- are analytic respectively in the upper and lower 
halves of the complex plane £ (similar features have been observed first in 
PH and later in [II [2]). 

The Jost eigenfunction (p and the constant eigenfunction p form a basis 
in the ring of solutions of Lf = 0; thus any solution / of Lf = is a 
function of tp and p. The analytic eigenfunctions if>±, in particular, possess 
the representations: 



the above analyticity properties of the LHS of (1251) in the complex £ - plane 
imply that x+(0 an d X-(0 are analytic respectively in the upper and lower 
halves of the complex plane £, decaying at £ ~ 00 like 0(£ _1 ). Therefore 
their Fourier transforms x+( u ) an d X-( u ) have support respectively on the 
positive and negative u> semi-axes. 

The spectral data x± can be constructed from the scattering datum a 
through the following linear integral equations 



ip± = £±(<p,p) = v + x±{<f,p), 



(24) 



defining the spectral data x±- 

Since the y — > — 00 limit of (12^1) reads: 

lim ip±-£= X±(t,P)> 



(25) 



X+(^,p) + 9(u) (a(u,p) + J R dr] x+(v,p)Q(v^,p)) = 0, 
X-(w,p) + B(-u)) (<j(w,p) + J R dr] x_(7],p)Q(r],uj,p)) = 0, 



(26) 



involving the Fourier transforms a and x± °f °" an d X± : 




(27) 



and the kernel: 




(28) 



To prove this result, one first evaluates (T24|) at y = +cxd, obtaining 




(29) 
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Applying the integral operator J R d^e~ lw ^- for uj > and uj < respectively 
to equations (l2T)j) + and (12T)1) _. using the above analyticity properties and the 
Fourier representations of x± an d cr, one obtains equations (1261) . 

The reality of the potential: i; G 38. implies that, for p £ ~R, Tp = tp, 
ip + = consequently: a = a, x+ = X— 



3.2 Inverse problem(s) 

As in [H12], due the ring property of the space of eigenf unctions, it is possible 
to construct three distinguished versions of the inverse problem, all based on 
equations (1241) . 

The first version is obtained subtracting £ from equations (1241) _ and ( 124"1) + . 
applying respectively the analyticity projectors P + and P_: 

P ± = ±— [ P—. (30) 

and adding up the resulting equations, to obtain the following nonlinear 
integral equation for the Jost eigenfunction ip: 

<p(x,y,p) + ^ j R p ,_%+ie) X- &(x,y,P'),p') - (31) 
A L P >-Z-ie) X+ (<P(?, y,p'),p') =x-py. 

Once ip is reconstructed from (13T1) . given x±(£,p), t ne analytic eigenfunctions 
follow from (T24|) . and the potential w from equation (122]) . This inversion 
procedure was first introduced in [18] and also used in [H [2] . 

The second version is the linear analogue of the nonlinear problem (13~Tj) . 
obtained exponentiating the Jost and analytic eigenfunctions used so far. 
Consider the following functions: 

^{x,y,p;a) = e ia ^ x ' y ' p \ tf ± (x, y,p; a) = e ia ^ x ' y ' p) , aGl (32) 

Due to the ring property of the space of eigenfunctions, also $(x, y,p; a) and 
ty±(x, y,p; a) are eigenfunctions; $(x, y,p;a) is a Jost eigenfunction charac- 
terized by the asymptotics $ — > exp(ia£), y — ► — oo, while \&±(x, a) 
are analytic respectively in the upper and lower halves of the p plane, with 
asymptotics: "$?± = exp(ia£)[l — (av)/p + 0{p~ 2 )]. 

Exponentiating the representations (T241) . one obtains the linear expan- 
sions of the analytic eigenfunctions ty± in terms of the Jost eigenfunction 
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$: 

%?±(x,y,p;a) = $(x,y,p;a) + f R d[3K±(a } (3 } p)$(x } y } p; (3), (33) 

where the new spectral data K± are defined in terms of \± °y : 

K ± (a,(3,p) = j R g e *(a-/»)€[ e fax±(e*) _ i]. (34) 

Multiplying the equations (l33l + and (l33l)_ by exp(— ia£,), subtracting 1, 
applying respectively P_ and P + , and adding the resulting equations, one 
obtains the following linear integral equation for $: 

$(x, y, p; a) + ± J R -^-^ J R d(3K_(a, (3, p')$(x, y, p'; (3)e^'- P )v- 



p' — (p+ie) 

f A 

p'-(p-ie) 



~hk W^hr) J R dpK + (a,P,p'Mx,y,p';P)e^'-^ = e*»C»-w). 



(35) 

Once $ is reconstructed from ( J35l) and, via (l33l) . \I/± are also known, the 
potentials are reconstructed in the usual way from the asymptotics of \I/±. 
The reality constraint (v = 0) implies, for a,p G R: = $(— a), 



*+(a) = Consequently: X+(a, /3) = X_(-a, 

The third version of the inverse problem is a Riemann problem. Solving 
the algebraic system (T24"|) _ with respect to <£> = (assuming local 

invertibility) and replacing this expression in the algebraic system (|24p + . one 
obtains the representation of the analytic eigenfunction ip + in terms of the 
analytic eigenfunction ?/>_: 

V>+ = W-,p)=V- + W-,p), peR. (36) 

where TZijp-.p) = tC + (C(i/;~,p),p), which defines a scalar nonlinear Riemann 
problem on the real p axis. The Riemann datum 1Z is therefore constructed 
from the data K,± by algebraic manipulation. Vice-versa, given R, one con- 
structs the solutions ip± of the nonlinear Riemann problem (13"oT) and, via the 
asymptotics (1221) . the potential t> . 

The reality constraint for 1Z takes the form: lZ(lZ(£,p),p) = £, V£, for 
pGK. 

Remark . Dressing schemes can be formulated from the three different inverse 
problems presented in this paper in a straightforward way. 
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3.3 t - evolution of the spectral data 

As v evolves in time according to ([TJ, the t-dependence of the spectral data 
a i X±i R an d K± is described by the following explicit formulas: 

a(£, p,t) = <r(Z- p% p, 0) , <S(£, P, t) = £ + <r(£ - p 2 t, p, 0) , 

x±(e,p, *) = x±(e - p%p, o), /c±(e,p, *) = e + x±(e - Ap, o), 

R(ti,p,t) = R(£ - p 2 t, p, 0) , p, t) = £ + - A P, 0) , ^ 

K±(a,(3,p,t) = K ± (a,/3,p,0)e^-^ 2 '. 

To prove it, we first observe that 

<t>(x,y,t,p) = ip(x,y } t } p) -pH (38) 

is a common Jost eigenfunction of L and M. The y = +oo limit of equation 
Mcf) = yields a t + P 2 &( = 0, whose solution is f[3"Th). Analogously, 

7T ± (.t, y , t, p) = ip ± (x,y,t,p) - pH (39) 
are common analytic eigenfunctions of L and M; therefore 

7r±=£±(0,p), 7T+ =7^(7T_,p), (40) 

for some functions KL± and 7?. depending on x,y,t only through and 7r_. 
Comparing at i = these equations with equations §M§ and (13T)|) . one ex- 
presses /C± and 7?. in terms of /C± and 7?., obtaining equations (13Tb. c). The 
t-evolution of .K± is obtained replacing fl37b) in 

Analoguosly, one can obtain the explicit t - dependence of the spectral 
data as v evolves according to the other equations of the hierarchy (J7J). 



4 One- parameter families of commuting 
dynamical systems 

It is well-known (see, f.L, [19]) that linear first order PDEs like ffT5j) . fl3|) . (fill 
are intimately related to systems of ordinary differential equations describing 
their characteristic curves. For instance, the dynamical systems associated 
with the one-parameter family of vector fields L, M in §S§ are: 

f=P + v x (x,y,t), 
M: ^=p 2 + pv x (x,y,t)-Vy{x,y,t). 
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Therefore equation (T7J) characterizes the class of functions v for which the 
two first order dynamical systems (JI ) commute Vp. 

There is also a deep connection between the above 1ST and the time 
(y) - scattering theory for the commuting flows (jUJ). Let (f>(x,y,t,p) be the 
common eigenfunc tion of L and M defined in (138]) ; then, solving the system 
uj = (f)(x,y,t,p) with respect to x (assuming local invertibility), one obtains 
the following common solution of PTT) : 

uj = (p(x,y,t;p) —p 2 t x = r(uj,y,t;p) ~ py+p 2 t + uj, y ~ — oo. (42) 

The y = +00 limit of the solution r(uj, y, t; p): 

x ~ py + p 2 t + Q(uj,p), y ~ +00 (43) 



defines the time (y) - scattering datum A(u,p) = £l(u,p) — uj of f|4Th). 
which is connected with the 1ST datum S by inverting the system uj = 
S(x — py — p 2 t,p, 0) with respect to x: 

uj = S(x — py — p 2 t,p,0) <^ x — py — p 2 t = Q(uj,p). (44) 

As a byproduct of the 1ST of this paper one can reconstruct, from the scatter- 
ing datum A(uj,p) of the one-parameter family of dynamical systems (f4Th ). 
the function v. 

Similar considerations can be made for the dynamical systems associated 
with the whole hierarchy of one-parameter family of operators L, M n , n G N. 
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